Local contact stress measurements at a rough interface 
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An original MEMS-based force sensing device is designed which allows to measure spatially re- 
solved normal and tangential stress fields at the base of an elastomeric film. This device is used 
for the study of the contact stress between a rough film and a smooth glass sphere under normal 
load. The measured profiles are compared to Finite Elements Method calculations for a smooth 
contact with boundary conditions obeying Amontons-Coulomb's friction law. The accuracy of the 
measurements allows to discriminate between dry and lubricated contact conditions and to evidence 
load-dependent deviations from Amontons-Coulomb's profiles. These deviations are qualitatively 
interpreted by taking into account the finite compliance of the micro-contact population. 

PACS numbers: 46.55.+d,81.40.Pq, 85.85.+J, 62.20.-x 



Knowledge of the stress field at the contact between 
two solids is of considerable interest in numerous con- 
texts such as mechanical engineering, seismology or solid 
friction. Starting with the classical analysis of Hertz this 
problem has been extensively studied theoretically in var- 
ious geometries for simple boundary conditions [1, 2, 3]. 
These calculations generally assume smooth interfaces 
which makes them unapplicable to most practical situ- 
ations: real solids are generally rough at the micrometer 
scale. The apparent contact is therefore made of a large 
number of isolated load-bearing micro-contacts, which 
is expected to modify the overall stress distribution as 
predicted in [4, 5]. Exhibiting such deviations requires 
to access the stress field with a spatial resolution be- 
low the contact size. In recent years, Micro Electro Me- 
chanical Systems (MEMS) have emerged as convenient 
tools for such stress measurements [6, 7]. Their quasi- 
monocristallinc Silicon structure offers a high linearity 
and a low hysteresis response. Their production pro- 
cesses, borrowed from micro-electronic technology, also 
allow for progressive miniaturization and production of 
arrays of sensors. 

In this Letter, we report on a new stress sensing device 
based on a triaxial MEMS force-sensor embedded at the 
base of rough, nominally flat elastomeric film (Fig. 1). 
The normal and tangential stress profiles are measured 
in the case of a spherc-on-plane contact under normal 
load. The high accuracy of the measurements allows for 
a direct comparison to elastic calculations under various 
mechanical boundary conditions. 

The elastomeric material is a cross-linked polydi- 
methylsiloxane (PDMS) (Sylgard 184, Dow Corning) of 
Young's modulus 2.2 ± QAM Pa, and of Poisson ra- 
tio 0.5 [8]. No measurable stress relaxation being ob- 
served after a sudden loading, the PDMS can be con- 
sidered as purely elastic. The film is obtained by pour- 
ing the cross-linker/PDMS melt on the MEMS into a 




FIG. 1: Schematics of the stress sensor. The MEMS sensi- 
tive part is a rigid cylinder (diameter 550 \im, length 475 /im) 
attached to a suspended circular Silicon membrane (radius 
1 mm, thickness 100 /im) whose deformation is measured via 
four couples of piezo-resistive gauges. It enables to measure 
the applied forces along three orthogonal directions at the 
base of a nominally flat rough PDMS film (thickness 2 mm, 
lateral dimensions 50 x 50 mm). The macroscopic normal and 
tangential loads applied on the film are measured through the 
extension of two orthogonal loading cantilevers (normal stiff- 
ness 641±5iV.m~ 1 , tangential stiffness 51100±700 N.m' 1 ) by 
capacitive position sensors (respectively MCC30 and MCC5, 
Fogale nanotech). 



parallelepipedic mold covered with a PolyMethylMeth- 
Acrylate plate roughened by abrasion with an aque- 
ous solution of Silicon Carbide powder (mean diame- 
ter 37 /urn). After curing and demolding the resulting 
rms surface roughness is measured with an interferen- 
tial optical profilometer (M3D, Fogale Nanotech) to be 
p = 1.82 ± 0.10 fim. This roughness is sufficient to avoid 
any measurable pull-off force against smooth glass sub- 
strates, as discussed in [9]. 

The force sensing device is first calibrated by indentat- 
ing the film surface with a rigid rod of diameter 500 fim 
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FIG. 2: Radial stress profiles under normal loading by a rigid 
thin rod. The measured tension U (•) is compared to FEM 
results (solid line) for a normal load of 1 N for (a) the normal 
and (b) the tangential stress. The error bars, when visible, 
represent the electronic noise. The black rectangular patch 
on (a) represents the rod diameter. 



(details of the complete calibration procedure are given 
in [10]). In this contact geometry the sensor measure- 
ment is found to be linear with the applied load. By 
successively varying the position of the contact with re- 
gards to the sensor symmetry axis along the y direction, 
one can reconstruct both the normal (zz subscript) and 
tangential (yz subscript) output voltage radial profiles U 
(Fig. 2). They are then compared to the results of Finite 
Elements Method (FEM, Software Castem 2007) calcu- 
lations for the stress a at the base of an axi-symmetrical 
elastic film (of same elastic moduli) perfectly adhering to 
its base and submitted to a prescribed normal displace- 
ment over a central circular area [19]. As expected for 
contacts of dimensions smaller than the film thickness, 
the stress calculated at the base of the film is found to be 
insensitive to the friction boundary conditions. For each 
stress component the relationship between the stress field 
a and the output tension U is assumed to be 



U(x,y) = A .G ® a(x,y) 



(1) 



where A is a conversion constant (expressed in mV/Pa), 
G is a normalized apparatus function and ® is a convo- 
lution product. 

G zz and G yz are obtained by numerical integration of 
Eq.l in Fourier space. The resulting functions exhibit 
a bell shape with typical width of the order of 600 /im 



comparable to the MEMS lateral dimension. G zz has 
proven to be adequate for all contact situations tested, 
therefore validating the hypothesis underlying Eq.l. In 
contrast the tangential output of the MEMS appears to 
be sensitive to the presence of normal stress gradients 
over the MEMS size, which prevents a solid determina- 
tion of G yz and hence robust measurements of the tan- 
gential stress [20]. Still, comparative information can 
be extracted when considering situations exhibiting low 
normal stress gradients, such as contacts against spheres 
having a large radius of curvature as studied in this Let- 
ter. For this subclass of situations an effective apparatus 
function having the same spatial distribution as G yz but 
an amplitude increased by 20% is used, for it provides sat- 
isfying agreement between experimental and calculated 
tangential stress profiles over the whole range of applied 
loads. 

Contact is made with an optical plano-convex spherical 
glass lens (radius of curvature 128.8 mm). Both the glass 
and the PDMS surfaces are passivated using a vapor- 
phase silanization procedure which reduces and homoge- 
nizes the surface energy. The contacts are obtained using 
the following loading sequence. The substrate is pressed 
against the film up to the prescribed load P within 2 % 
relative error. Due to the associated tangential displace- 
ment of the extremity of the normal cantilever, a sig- 
nificant tangential load Q is induced. From this posi- 
tion, the contact is renewed by manual separation which 
results in a much smaller but finite Q. The substrate 
is eventually translated a few \im tangentially down to 
Q = 0. Both the surface treatment and the loading se- 
quence yield an excellent measurement reproducibility. 
As for the rod indentation the radial profiles are derived 
assuming a homogeneity of the surface properties of the 
film from a series of 33 contacts whose centers lie every 
0.5 mm along the y direction. These profiles divided by 
Aq are denoted as S and have the dimension of a stress. 

For quantitative comparisons, FEM calculations are 
carried out for a frictional sphere-on-plane contact with 
the same geometry. Both contacting surfaces are taken 
to be smooth. The interface is assumed to obey a local 
Amontons-Coulomb's law of friction which states that no 
slip occurs wherever the ratio of tangential over normal 
stress remains smaller than the static friction coefficient 
IX. The normal displacement of the rigid elastic sphere 
is prescribed. Both solids are discretized with a uniform 
mesh size of 50 \im and the normal displacement of the 
rigid elastic sphere is prescribed. The contact conditions 
are satisfied using a double Lagrange multiplier implying 
that both surfaces are slave and master. The normal load 
is reached step-by-step: at each step an iterative Newton- 
Raphson method is used to satisfy both the unilateral 
contact and the friction law. 

Figure 3 compares the S and G <g) <J profiles for three 
values of P. Within the error bar the measured profiles 
are bounded by both numerical profiles obtained with 
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FIG. 3: Measured stress profiles (S) under normal loading 
by a rigid sphere (P = 0.34 N (•), 1.37 N (■), 2.75 N (a)). 
Comparison is made with G ® a for fi — (solid line) and 
fi = oo (dashed line), (a), (b) Normal stress in linear and 
semi-logarithmic scale respectively, (c) Tangential stress. The 
black rectangular patches on (a) represent the contact diam- 
eters obtained from the FEM calculations for fi = 0. 



fi = and n = oo over the whole spatial range and over 
3 orders of magnitude, as clearly displayed on Fig. 3b. In 
the outer region of the contact S is systematically very 
close to the frictionless profile whereas at the center it 
increasingly departs from it with the load, as discussed 
further. The slight asymmetry of the tangential stress 
profiles (Fig. 3c) can be related to some irreversible micro- 
slip which occurs during the loading sequence. 

The effect of friction can be probed by performing sim- 
ilar measurements when the contact is lubricated. A 
glycerol droplet is inserted at the interface prior to load- 
ing, allowing for a complete relaxation of the interfacial 
tangential stress. This is done for two limiting loads 
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FIG. 4: Measured normal stress profile (S) under normal load- 
ing (P — 2.75 N) by a rigid sphere for both dry (•) glycerol 
lubricated (■) contacts. Shown in solid and dashed lines are 
the G <8> a profiles for 4 values of fj, (0, 0.3, 0.4, oo). The 
black rectangular patch represent the contact diameter ob- 
tained from the FEM calculations for fj, = 0. 
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FIG. 5: Load dependence of the maximum normal stress Smax 
at the base of the elastic film as a function of P, for dry 
(•) and lubricated (■) contacts), together with the maximum 
calculated normal stress Pmax for various values of fi. All 
values are normalized by p^ M . 



(P = 0.69 AT and P = 2.75 N). The lubricated profiles 
are characterized by both a larger spatial extent and a 
lower maximum amplitude than what is found for the 
dry ones. In addition the lubricated profiles are found to 
follow almost perfectly the frictionless G <S> a curves (as 
shown for S zz on Fig. 4 at P = 2.75 N). 

For dry contacts and sufficiently high loads, we find 
similarly that the normal stress profiles can be strictly 
bounded by the calculated ones for two very close values 
of fx, hence providing an interval for an effective friction 
coefficient fi e (for instance for P = 2.75 N, fj, e is in the 
[0.3, 0.4] range as shown on Fig. 4) which value is found 
to increase with P. However, for small enough loads 
the experimental profiles are beyond the possible range 
for smooth contacts and [i e cannot be defined. This be- 
havior can be well caught by plotting for each load the 
maximum value of the pressure profiles S ma x normalized 
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by Pmax = G Z z <£> cr(0) expected for a smooth frictionlcss 
contact (Fig. 5). It is compared to P^ax/Pmax calculated 
for fi varying from to oo. Whenever possible the value 
of /i e can be directly read from Fig. 5. Clearly the exper- 
imental data cannot be described within the Amontons- 
Coulomb's framework using a single friction coefficient. 

The values of p e can be compared to the actual static 
friction coefficient [i s defined as the ratio of the tangential 
over the normal loads at the onset of macroscopic sliding. 
fi s has been measured for a driving velocity v = lOOfim/s 
of the rigid base of the sensor and is found to decrease 
from 1.8 to 1.5 as the load is increased from 0.34 A to 
2.75 A, a behavior which is usually attributed to the fi- 
nite adhesion energy of the interface [11]. We observe an 
opposite load-dependence for /i e , whose understanding 
implies therefore to reconsider the Amontons-Coulomb's 
assumption of a smooth undeformable interface. 

A description of the normal compressibility of a rough 
interface [12, 13, 14] has first been proposed by Green- 
wood and Williamson [15] by modeling the interface as 
an ensemble of spherical asperities of equal radius whose 
heights are statistically distributed. They showed that 
the true area of contact is proportional to the confin- 
ing load thus providing an understanding of the quasi- 
indcpcndence of the static friction coefficient with the 
apparent area of contact. When adapted to a sphere-on- 
plane contact [4] this model predicts an increase of the 
apparent radius of contact with respect to Hertz's theory 
and a decrease in the maximum normal stress. These 
effects are expected to be enhanced for small loads as 
the normal displacement of the rigid substrate becomes 
comparable to the thickness p of the rough layer. This 
could explain why S ma x/Pmax ^ aus bcl° w 1 for vanish- 
ingly small applied loads. 

At high enough load, the fact that the lubricated pro- 
files are very close to those calculated for a smooth fric- 
tionlcss contact suggests that these normal compressibil- 
ity effects are negligible. The finite shear stiffness of the 
interface, evidenced experimentally in [16, 17], is then ex- 
pected to be primarily responsible for the drop from p s 
to fi e . It allows for a partial relaxation of the tangential 
stress before the onset of slippage, which qualitatively 
translates into an apparent reduction of the friction co- 
efficient. This reduction should be larger the smaller the 
load, when the relaxed tangential displacement becomes 
comparable to the slip amplitude expected for a friction- 
less smooth contact. This behavior is consistent with 
what is observed on Fig. 5. At this point a quantitative 
description of the mechanical state of the system should 
involve an elasto-plastic like friction law accounting for 
the ability of the asperities to deform reversibly before 
slippage. 

The stress sensor described in this Letter has proven to 
be well-suited for the study of stress fields in centimeter- 
sized contacts. The measurements are accurate enough 



to discriminate between lubricated and dry contacts and 
to evidence deviations from Amontons-Coulomb's model 
of friction induced by a micron-sized interfacial rough- 
ness. These measurements can then provide a fine test 
for any mechanical model for the frictional interface. 
Many other aspects of friction and contact mechanics, 
such as the dynamical frictional regimes or the history- 
dependence of a contact submitted to an oscillatory tan- 
gential load below the sliding threshold, are likely to be 
probed with the same device by applying controlled tan- 
gential load sequences. This device might also be of in- 
terest in domains such as rhcology or adhesion where 
accurate spatially resolved stress measurements at inter- 
faces are needed. 
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discussions. 
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